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AMENDMENTS TO THE C L AIMS : 

This listing of claims will replace all prior 
versions, and listings of claims in the application: 

L I STING OF CLAIMS : 

1. (cancelled) 

2. (currently amended) A method according to claim 
8—1-, for simulating an oscillating phenomenon, wherein 
the digital model comprises, on the one hand at least one 
linear part (3), representing the input impedance or 
admittance of the resonator, and, on the other hand, one 
non-linear part (2) modelization the role of the 
excitation source (22) of the phenomenon to be simulated. 

3. (currently amended) A simulation method according to 
claim 8— 3r, for real time digital synthesis of an oscillating 
phenomenon, wherein, from a system of equations between at least 
two variables representative of the behaviour of a wave in the 
resonator, an expression of the input impedance or admittance of 
the resonator is established in the form of a linear filter 
including delays, without any decomposition into two-way waves, in 
order to realise at least one linear part (3) of the model. 

4. (original) A method according to claim 3, wherein the 
linear part (3) of the model is coupled with one non-linear part 
(2) involving the evolution of the non-linearity as expressed 
between the two variables of the input impedance or admittance 
relation of the resonator. 

5. (original) A method according to claim 4, wherein the 
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linear part (3) of the digital simulation model of the impedance 
or admittance equation in based on to two elementary waveguides 
fulfilling a transfer function between the two variables of the 
impedance or admittance relation. 

6. (original) A method according to claim 5, wherein the 
linear part (3) with two waveguides of the model is coupled with a 
loop connecting the output to the input of said linear part (3) 
and comprising a function (21) involving the non-linearity as 
expressed physically. 

7. (original) A method according to claim 6, wherein the 
model is driven by at least two parameters representative of the 
non-linear physical interaction between the source and the 
resonator, by means of a loop connecting the output to the input 
of the linear part (3) and comprising a non-linear function (21) 
playing the part of an excitation source for the resonator. 

8. (currently amended) a digital simulation method of a 
non : line ar interaction betwe en an excitation source and_^_wave_in 
a resonat or, by means of digital signal calculation toqls_based_on 
eguatign 3_the solution of which corresp o nds to the physical event 
of g_ phenome non to be simulated which can t ra nslat ed £ _at_each time 
and at each poin^_gX^the_J^3gnato^ 

two variables representative of the ef fect and of thg_„gguse^f 
said phenomenon to be simulated comprising: 

t-r anscribinq the impedance or admittance equation directly 
into a .digita l model form enabling to realise a non- linear 
intera ction between the two variab les_of__th e impedance ,jr 

admittance relation, 

A m etho d i ' ' " J i d nim 1 wherein said method ^s^dapted 

for real-time sound synthesis of a musical instrument comprising, 
at least, one excitation source with non-linear characteristics 
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and a linear resonator, the sound produced by the instrument 
resulting from a coupling, between the excitation source and the 
resonator, expressed at least by a linear impedance or admittance 
relation and a non-linear relation between two physical variables 
representative of the effect and of the cause of the sound 
produced, a method where the sound produced by the instrument is 
simulated, in real time, by modeiization the physical phenomena 
governing the operation of the instrument, wherein, to realise 
said physical modelization, the impedance or admittance linear 
relation is expressed directly and digitally between two physical 
variables representative of the cause and of the effect of the 
phenomenon to be simulated and said impedance or admittance 
relation in digital form is associated with the non-linear 
relation between the same variables. 

9. (original) A method according to claim 8 for synthesis of 
the sound of an instrument with complex resonator, wherein the 
resonator is decomposed into a series of successive elements and 
wherein the impedance or admittance relations corresponding 
respectively to each element of the resonator are calculated and 
combined xn order to obtain a global impedance corresponding to 
the geometry of the resonator. 

10. (previously presented) A method according to claim 8, 
wherein, for real-time synthesis of the sound produced by a wind 
instrument, the two variables of the impedance relation are the 
acoustic pressure (p e ) and flow (u e ) at the input of the resonator. 

11. (original) A method according to claim 10, wherein, for 
an open-ended cylindrical resonator, the linear part (3) of the 
digital transcription model of the impedance equation is the sum 
of two elementary waveguides having as excitation source the flow 
(u e ) at the input of the resonator, and fulfils the transfer 
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function 



7 Pe(Q3) = 1 exp(-~2ik(o)L) 

e ^°' Ue(co) 1 + exp(-2ik(a>)L) 1 + exp(-2ik(co)L) 



(12) 



wherein : 

co is the angular frequency of the wave ; 

Ze(to) is the input impedance of the resonator, 

Pe(co) and Ue(co) are the Fourier transforms of the 
dirnens i onles s values of the pressure and of the flow at 
the input of the resonator ; 

k(<30) is a function of the angular frequency which 
depends on the phenomenon to be simulated ; 

L is the length of the resonator. 

12. (original) A method according to claim ll r wherein each 
of the two waveguides involves a filter having as a transfer 
function : 

- F{g>) 2 - - exp(-2ik{o)L) 

and representing a two-way travel of a wave, with a sign 
change at the open end of the resonator, each waveguide 
corresponding to a term of the impedance equation. 

13. (original) A method according to claim 12, wherein the 
model is driven by the length ( L) of the resonator and at least 
two parameters (£, y) representative of the non-linear physical 
interaction between the pressure (p e ) and the flow (u e ) at the 
input of the resonator, by means of a loop connecting the output 
to the input of the linear part (3) and comprising a non-linear 
function (21) as an excitation source for the resonator. 
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14. (original) A method according to claim 11, wherein the 
non-linear function has the pressure and the displacement of the 
vibration formation member as input parameters, and is controlled 
by at least two parameters simulating a player playing. 

15. (original) A method according to claim 14, 
wherein the playing parameters for controlling the non- 
linear function are: 

a parameter q characteristic of the mouthpiece and of 
the action of the player on the vibration formation 
member , 

a parameter y representative of the pressure applied 
to the vibration formation member. 

16. (previously presented) A method according to claim 10, 
wherein, for real-time synthesis of the sound to be simulated, a 
formulation is realised in the time domain of the angular 
frequency response of the impedance of the resonator, by 
approximation of the losses represented by the filter by means of 
an approximated digital filter. 

17. (original) A method according to claim 16, wherein, to 
express the angular frequency response of the impedance of the 
resonator, a one pole digital filter is used of the form: 

Si ^ b 0 exp(-2teD) 
r{w) = - — — r 

1 - a 1 exp(Htn) 

(13) 

wherein : 

CO 

T e f / e being the sampling frequency, 
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D = /. - 

c is the pure delay corresponding to an away or 

return travel of the wave in the resonator, 

the coefficients bo and al are expressed in 

I |2 |~ ]2 

relation to the physical parameters so that ]F(ro) 2 l = [F(to)[ 
for a value 

col of the angular frequency corresponding to 
the fundamental playing frequency and another value co2 
corresponding to a harmonic, 

and the following differential equation is derived: 

p e (n) - u e (n) - a,u e (n - 1) - b D u e (n - 2D) + a.pjn - 1) - b oPe (n - 2D) (16) 

18. (original) A method according to claim 17, 
wherein the coefficients bo and al are obtained by 
solving the equation system: 

|FK) 2 | 2 (1 + a 2 - 2a, cos(tnj) - b 2 
|F((o 2 ) 2 j 2 (1 + a 2 - 2a, cos(m 2 )) - b 2 0 

with |F(03) 2 | 2 = exp(-2aC^L), said coefficients being given 
by the formulae: 



a, - 



A1 - A2 - J(A71 A 2 ) 2 -(F, -F 2 ) 2 



2 



^F.F^c, - c^A, - A 2 - V(A, - A 2 ) 2 - (F, - F 2 ) r ) 



j o F - F 



2 
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wherein 

c 1 = cosCtu^. c 2 = 005(1352), F t = jF(ot) 1 ) 2 | 2 ,F 2 = |F(co 2 ) 2 j 2 J 
A, - FjC.,, A 2 = F 2 c 2 

19. (original) A method according to claim 18, for 
simulating a cylindrical resonator instrument, from a physical 
modelization governed by the system of equations : 

1 ^) + %^) + x(t) = ± p e (t) 
a) 2 dt 2 co r dt 

(with the sign + for a reed and the sign - for the 

lips) 



P.(a>) = i tanf— - ~ ac» 1/2 L\j e (co) 

u e (t) = 1(1 - sign( Y - x(t) - 1))sign(y - p a (t)K(1 - y + x(t))^/|y - p.(t)| 



wherei 



n cor is the resonance frequency and qr is the quality 
factor of the reed or of the lips, wherein that said 
system of equations is solved in the time domain from an 
equivalent sampled formulation of the angular frequency 
response of the displacement of the reed or of the lips 
and of the impedance relation which is translated by the 
system of equations : 

x (n) =bi a pe (n-1) +a Xa x (n-1) +a 2a x (n-2 ) 

(18) 

p e (n) =u e (n) -aiu 0 (n-1) -b 0 u a (n-2D) +aip e (n-1 ) b 0 p s (n~2D) 
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A 

u e (n) - - (1 - sign(y - x(n) - 1))sign(y - p.<n))?n)- y + x(n))^/[Y - p e (n)| (20) 

said equations being used sequentially by grouping the terms 
not depending on the time sample n, in order to calculate in 
succession : 

x (n) =b ia p e (n-1) +a ia x (n-1) +a 2a x (n-2) 

(21) 

V— aiu e (n-1) -b 0 u e (n-2D) +ajp e (n-1) -b 0 p e (n-2D) 

(22) 

W = -1(1 - sign( Y - x(n) - 1))Q) - Y + x(n)) 

(23) 

u e (n) - ^sign(Y - V)(-bc 0 W 2 + W^/(bc 0 W) 2 + 4|y -~V| (24) 
p G (n) -b 0 c a u e (n) +V 

(25) 

20. (original) A method according to claim 19, for more 
realistic simulation of the produced sound, wherein, by neglecting 
the radiation, the external pressure is expressed as the time 
derivation of the outgoing flow, in the form: 

P«i(t) = ^(p.(t) + u a (t)) 

(26) 

and is calculated, at each sampled time (n) , by differential 
between the sums of the internal pressure pe and of the flow ue, 
respectively at the time (n) and at time (n-1) „ 
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21. (original) A method according to claim 19 for simulating 
a multimode reed instrument, wherein the calculation of the 
acoustic pressure and of the flow at the mouthpiece is performed 
by sequential resolution of a system of equations wherein the 
displacement of the reed at each time(n) is in the form: 

x (n) ^b a ip e (n-l) +b a2 p e (n-2) +b aD1 p e (n-D a -l) +a ai x (n-1) n-a a2 x (n- 
2 ) +a aD x (n-D a ) +a. aDi x (n-D a ™i ) 

the coefficients aal, aa2, aaD2, aaDl being defined 

by: 



and the coefficients bal, ba2 , baDl by : 



C - Ca a . - Cb, 



fe fe fe 



(28) 

2 1 

by setting : A 1 - . = , A 2 -;-co r q r and A 3 =A 2 A L q r , 

co rA /q r 2 +4 * 

the following equations being the same as for a 

single mode reed. 



22. (previously presented) A method according to 
one claim 19, wherein that a model is prepared for a 
cylindrical resonator with terminal impedance from the 
basic model corresponding to a cylindrical resonator and 
being the sum of two waveguides involving each a filter 
having as a transfer function -F(w) 2 --exp{~2ik(co)L) , 
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while replacing the expression exp ( -2 ik ( (D ) L) by the 
expre s s ion 

D/ x Z c -Z s (co) 

R (co) exp (-2ik (©) L) , wherein R(ra) = - 

Zc being the characteristic impedance 

DC 

p etZ s the output impedance 



nR 2 5 " U B (<o) 

23. (previously presented) A method according to claim 19, 
wherein, from the impedance model for cylindrical resonator 
instrument and the associated differential equations, other more 
complex impedance models are built for simulating oscillating 
phenomena produced by a resonator of any shape by combining 
impedance elements in parallel or in series and by using digital 
approximations for an explicit use of the physical variables 
involved in the production of said oscillating phenomena and a 
more flexible control of the result of the simulation. 

24, (previously presented) A method according to claim 23, 
wherein that, from the basic model of a cylindrical resonator 
wherein the angular frequency response of the displacement of the 
reed or of the lips is translated by a system of differential 
equations providing, at each time (n) , the displacement x(n) the 
pressure pe(n) and the flow ue{n) at the input of the resonator, a 
model for a conical resonator is built wherein the equation of the 
pressure is in the form: 

p e ( n ) - bc 0 u e ( n ) +bciue ( n — 1 ) +bc 2 u e ( n-2 ) +bc D u e { n-2 D ) +bc D iU e ( n- 2 D- 

1) 

+ acipe (n-1 ) +ac 2 pe (n-2) +ac D p e (n-2D) +ac D ip e (n-2D-l) 

(33) 

wherein the coefficients bcO, bcl, bc2 , bcD and bcDl 
are defined by : 
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bCr 



a 1 + 1 



bc 2 = 



foc D = 



V 01 



OC D i=- 



and the coefficients aci, ac 2 , ac D and ac Di are defined by : 



ac 1 = 



a,G p + G r 



ac 2 = 



ac D = 



bo G m 



ac D1 - b 0 



by noting : G p - 1+ 



1 



2f e ^ 
e c 



and G„ 



2f ^ 



25. (previously presented) A method according to claim 2 4 
wherein that, from the basic model for a cylindrical resonator, a 
model for a short resonator having a length 1 is built, by an 
approximation of the impedance according to the expression: 



tan (k{co)l) = G(co) + icQH(co) 



(34) 



w 



herein G(co) = 



1 - exp - act. 

CO 



1 + exp) - ac^f - I 



and H (co) =-(1 - G(a>)) . 
c 



26. (previously presented) A method according to 
claim 25 for simulating a wind instrument, wherein the 
mouthpiece or the bill is modelled by a Helmholtz 
resonator comprising a hemispheric cavity coupled with a 
short cylindrical pipe and a main resonator with a 
conical pipe, the input impedance of the resonator 
assembly which may be expressed as: 



13 



Docket No. 0552-1014 
Appln. No. 10/533,336 



Z e (<D) = ^ 

ico h 



1 

z' 



P° 2 io>^ 



ico — — 



1-f 



itan(k 2 (co)L 2 ) 



wherein V=^dR^ is the volume of the hemispheric 
6 

cavity, Ll is the length of the short pipe, L2 is the 
length of the conical pipe, Zl and 22 are the 
characteristic impedances of both pipes which depend on 
their radii, kl(co) and Jc2(a>) take into account the 
losses and the radius Rl and R2 of each pipe, and that, 
from the basic model for cylindrical resonator, from its 
extensions to the conical pipe and to the short pipe, a 
resonator model is prepared by expressing the pressure 
at the mouthpiece or in the bill by the differential 
equation : 

Pe (n) = ]Tbc k u e (n - k) + "f bc Dk u e (n - k - 2D) (36) 
+ k |>c K p e (n-k) + K f>c Dk p e (n~k -2D) 

k=1 *=o 

27. (currently amended) A method according to claim__8— 1, for 
the simulation of an oscillating phenomenon wherein both physical 
variables of the linear relation are the strength applied to one 
point of a mechanical system such as a string generating 
vibrations and the speed at this point, wherein that the 
admittance is expressed, at this point, in the form of a 
combination of the admittances of each part of the string, on both 
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sides of said point, each mechanical admittance being obtained 
from the basic model describing the acoustic impedance of a 
cylindrical pipe resonator, by expressing the speed at the point 
considered of the string in relation to the strength applied to 
this point, where the filter F(co) of the basic model may be 
expressed from a bending wave propagation model in a having a 
stiff eness . 

28. (previously presented) A digital device to implement the 
method according to claim 8, for the simulation of a musical 
instrument generating a sound resulting from a coupling between a 
Unear resonator and an excitation source with a non-linear 
characteristic, which can be expressed at least by a linear 
impedance or admittance relation and a non-linear relation between 
two variables representative of the effect and of the cause of the 
produced sound, said simulation device comprising a control 
element 1 including at least one gestural sensor 1 transforming 
the actions of a player into control parameters, a modelization 
element II including one non-linear part (2) associated with a 
linear part (3) and an element creating the synthesised sound III, 
wherein that the linear part (3) includes a computing block (31) 
driven by the length (L) of the resonator, having as an input 
parameter a signal representative of one of the variables, cause 
or effect, computed by the non-linear part (2), and the transfer 
function of which is the input impedance or admittance of the 
resonator, in that the non-linear part (2) implements a non-linear 
function (21) driven by at least two control parameters and having 
as input parameters a signal representative of the other variable, 
cause or effect, computed by the linear part (3) and a signal 
modelization the role of the excitation source, the linear part 
(3) being thus coupled with the non-linear part (2) in a closed 
loop, and in that the element creating the sound III computes a 
sound signal from signals representative of the cause and of the 
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effect of the sound to foe simulated, emitted respectively by the 
linear part (3) and the non-linear part (2) , 

29. (previously presented) A method according to claim 19, 
wherein that, from the basic model of a cylindrical resonator 
wherein the angular frequency response of the displacement of the 
reed or of the lips is translated by a system of differential 
equations providing, at each time (n) , the displacement x(n) the 
pressure pe (n) and the flow ue (n) at the input of the resonator, a 
model for a conical resonator is built wherein the equation of the 
pressure is in the form: 

p e (n) - bc 0 u e (n} +bc x u e (n-1) +bc 2 u e (n-2 ) +bc D u e (n^2D) +bc D iu e (n-2D- 

1) 

+ ac^e (n-1) +ac 2 Pe (n-2) +ac D p e (n~2D) +ac D ip e (n-2D-l) 

(33) 

wherein the coefficients bcO, bcl, bc2, bcD and bcDl 
are defined by : 

bc 0 = , bc-r - -~ f bc 2 =-^, dc d = - — - , bc D i- — 

G p G p G p Cs p u p 

and the coefficients aci, ac 2 , ac D and ac Di are defined by : 



p 



1 , n 1 

by noting : G D = 1+ — ana b ra - I 

2f. ** 2f 0 c 

6 c c 



30. (previously presented) A method according to claim 19 
wherein that, from the basic model for a cylindrical resonator, a 
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model for a short resonator having a length 1 is built, by an 
approximation of the impedance according to the expression: 



Zi(<b)= i tan (k«o)l) == G (co) + icoH(co) 



(34) 



1 - exp 



wherein G(co) = 



1 + expl - ac^| — I 



and H (co) - G(co)) 

c 
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